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ABSTRACT 

A  Clough-Tocher  like  interpolation  scheme  is  described  for  values  of 
position,  gradient  and  Hessian  at  scattered  points  in  two  variables. 

The  domain  is  assumed  to  have  been  triangulated.  The  interpolant  has 
local  support,  is  globally  twice  differentiable,  piecewise  polynomial , 
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SIGNIFICANCE  AND  EXPLANATION 


So—  design  applications  (e.g.  in  the  aircraft  industry)  require 
twice  differentiable  surfaces.  In  this  paper,  an  interpolation  scheme 
for  the  construction  of  such  surfaces  is  derived. 

More  specifically,  the  schene  requires  values  of  position,  gradient, 
and  Hessian  at  scattered  points  in  two  variables.  The  domain  is  assumed 
to  have  been  triangulated.  The  Interpol ant  has  local  support  (i.e. 
evaluation  at  a  point  in  a  specific  triangle  requires  data  only  on  that 
triangle),  is  globally  twice  differentiable,  piecewise  polynomial,  and 
reproduces  polynomials  of  degree  up  to  three  exactly.  Explicit  formulas 
are  given  for  the  coefficients  of  the  interpolant.  A  pilot  (FORTRAN) 
code  is  available  from  the  author. 


ilbility  for  the  wording  and  views  expressed  in  this  descriptive 
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A  BIVARIATE  C2  CLOUGH-TOCHER  SCHEME 
Peter  Alfeld* 

1.  Introduction 

In  order  to  save  space,  the  reader  is  assumed  to  be  familiar  with  the 

introduction  and  section  2.  of  Alfeld,  1984.  There,  a  bivariate  and  a 

trivariate  C1  Clough-Tocher  scheme  were  derived.  Having  obtained  a  C1 

scheme,  it  is  natural  to  use  the  same  techniques  to  construct  a  (bivariate) 

C2  scheme.  Such  a  scheme  would  be  useful  e.g.  for  the  design  of  aircrafts, 

2 

where  C  smoothness  of  surfaces  is  required.  Interpolation  by  a  function 

that  is  polynosdal  on  each  macrotriangle  would  require  a  polynomial  of  degree 

at  least  nine  and  data  through  fourth  order  (Zenisek,  1970).  In  this  paper, 

we  derive  a  scheme  that  is  piecewise  quintic  on  each  macrotriangle  and  that 
2 

requires  only  C  data. 

Parin,  1980,  has  shown  that  a  piecewise  polynomial  C2  scheme  cannot  be 

constructed  on  the  simple  Clough-Tocher  split.  The  difficulty  can  be  traced 

to  the  fact  that  each  vertex  angle  is  divided  into  only  two  parts.  For  a 
2 

C  scheme,  a  division  into  at  least  three  subangles  is  required.  The  work 

described  in  this  paper  grew  out  of  efforts  by  Arner,  Barnhill,  Farin,  and 

2 

Little  to  construct  a  C  scheme  for  the  split  depicted  in  figure  1.  That 
split  generates  the  minimum  number  of  triangles  (seven)  while  trisecting  each 
angle  at  the  vertices.  It  is  still  an  unsettled  question  if  a  C  scheme 
exists  for  the  minimal  split. 
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Fi*ure  1:  A  Minimal  SpUt 


In  the  approach  described  here,  the  Clough-Tocher  split  is  applied 
twice.  First  we  divide  the  macrotriangle  into  three  subtriangles ,  and  then  we 
divide  each  subtriangle  into  three  microtriangles .  Thus  each  macrotriangle  is 


the  union  of  nine  microtriangles.  We  refer  to  the  centroid  of  the 

macrotriangle  simply  as  the  centroid,  and  to  the  centroids  of  the  subtriangles 

as  the  subcentroids.  The  construction  is  illustrated  in  figure  2.  We  assume 

we  are  given  C  data  (i.e.  values  of  position,  gradient  and  Hessian)  at  the 

vertices,  and  construct  an  interpolant  that  is  quintic  on  each  microtriangle 

and  twice  differentiable  everywhere  on  the  triangulated  domain. 

The  major  difference  between  the  approach  in  Alfeld,  1984,  and  that 

described  here  is  that  for  the  C1  schemes  it  was  possible  to  identify  in  a 

natural  way  one  Bezier  ordinate  with  each  of  the  conditions  defining  the 

2 

interpolant.  This  turns  out  to  be  impossible  in  the  C  case. 
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2.0  Notation  and  Internal  Continuity 


We  consider  a  general  macrotriangle  with  vertices  V2,  and  V3, 

and  denote  the  centroid  by  V4  *  (v1  +  V2  +  V3)/3.  The  subcentroids 
are  denoted  by  V5  *  (V2  +  V3  +  V4)/3,  Vfi  *  (V,  +  v3  +  V4)/3,  and 
V7  ”  (V,  +  V2  +  V4)/3.  Figure  2  illustrates  the  construction. 

It  is  significant  that  the  triples  vi*v4*v5  •  v2,v4'v6  '  and 

V3'V4'V7  are  each  colinear.  This  can  be  seen  for  example  by  writing 

V5  -  <V4  +  V2  +  V3)/3 

-  (V4  +  (3V4  -  V1  -  V3)  +  (3V4  -  V,  -  V2))/3 

-  (7V4  “  (V,  +  V2  +  V3)  -  V, )/3 

-  (4V4  -  V, )/3. 

The  technique  of  expressing  the  location  of  a  point  in  terms  of  different 
sets  of  internal  and  external  vertices  will  be  used  frequently  in  the 
sequel.  To  obtain  smoothness  conditions  we  will  express  cross-boundary 
directions  in  terms  of  the  vertices  of  each  of  the  adjacent  microtriangles, 
thereby  facilitating  differentiation  on  each  of  the  microtriangles.  This 
approach  is  described  in  detail  in  Alfeld,  1984,  and  is  further  exenq>lified  in 
section  2.2.2  below. 

The  location  of  a  general  point  P  is  expressed  as 

7 

P  »  l  b  V 
i-1  1  1 

where  the  b^  are  the  piecewise  linear  cardinal  functions  defined  on  the 
triangulation  of  the  macrotriangle  by  b^(V^)  *  6^,  ®  being  the  Kronecker 

Delta.  The  interpolant  to  be  constructed  is  of  the  form 


51 


p(P) 


7 

n  (i.t) 
3-i  3 


7 

n 

3-1 


where,  by  convention,  0®  :=  1. 

The  function  p  is  continuous  on  the  macrotriangle  and  contains  121  free 
parameters.  (Remarkably,  although  to  all  appearances  coincidentally,  this  is 
the  same  number  of  parameters  as  is  available  for  the  construction  of  the 
trivariate  C1  Clough-Tocher  scheme,  see  Alfeld,  1984.)  For  a  given  point 
P,  the  generalized  barycentric  coordinates  can  be  easily  computed  as  follows: 

1.  Compute  the  barycentric  coordinates  with  respect  to  the 

macrotriangle,  i.e.  write  P  *  +  B2V2  +  ®3V3*  Let  {i,j,k}  =  {1,2,3} 

and  suppose  <  Bj  and  B^  <  B^.  Then  P  lies  in  the  subtriangle  with 
vertices  Vj ,  ,  and  V^.  The  barycentric  coordinates  with  respect  to  that 

triangle,  aj,  afc,  and  a4,  say,  are  given  by  a4  **  3Bif  aj  *  Bj  -  B±, 

ak  *  ®k  "  Bi* 

2.  After  relabeling,  repeat  step  1  on  the  subtriangle  to  locate  the 


microtriangle  containing  P. 


The  microtriangle  with  vertices  Vi#  Vj,  and  will  be  denoted  by 


T 


ijk* 


2. 1  Interpolation  to  Vertex  Data 

Since  it  is  impossible  to  associate  in  a  natural  and  unique  way  Bezier 
ordinates  with  conditions  on  the  interpolant  we  start  by  considering 
interpolation  to  the  vertex  data.  This  will  keep  a  larger  part  of  the  linear 
system  triangular.  Thus  we  do  not  assume  that  we  have  already  enforced 

2 

smoothness.  On  the  contrary,  interpolation  at  the  vertices  will  force  C 
smoothness  at  the  vertices.  It  affects  all  Bezier  ordinates  that  are  no  more 
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than  two  levels  removed  from  the  vertices.  A  simple  count  using  figure  2  or  3 
shows  that  45  conditions  have  to  be  satisfied.  These  are  given  as  equations 
1-45  in  the  appendix.  The  technique  for  the  derivation  is  analogous  to  that 
in  Alfeld,  1984,  and  in  Barnhill  and  Farin,  1981.  Note  that  the  data  are 
expressed  as  derivatives  in  the  direction  of  edges.  They  have  to  be  supplied 
consistently.  This  is  automatic  if  all  such  data  are  computed  f*'  m  cartesian 
derivatives  through  second  order. 

2.2  Internal  First  Order  Differentiability 

To  force  differentiability  of  the  interpolant  everywhere  in  the 
macrotriangle,  we  proceed  in  stages,  distinguishing  between  differentiability 
in  the  interior  of  subtriangles,  differentiability  across  lines  separating 
subtriangles,  and  differentiability  at  the  centroid.  The  complete  set  of 
conditions  is  described  in  figure  3.  Each  stencil  gives  coefficients  of  a 
linear  combination  of  the  Bezier  ordinates  that  must  be  zero.  For  each  type, 
one  representative  stencil  is  depicted.  For  better  legibility,  it  has  been 
moved  out  of  the  drawing  of  the  macrotriangle.  All  other  stencils  of  the  same 
type  can  be  obtained  by  drawing  them  about  the  appropriately  marked  centers. 
Stencils  in  subsequent  figures  should  be  interpreted  similarly. 

2.2.1  C1  Conditions  in  the  Interior  of  Subtriangles 

As  was  shown  in  Alfeld,  1984,  the  relevant  condition  is  that  the  Bezier 
ordinates  in  the  interior  of  internal  edges  equal  the  averages  of  the  three 
neighboring  vertices.  However,  interpolation  at  the  vertices  already  imposed 
two  conditions  on  each  internal  edge  emanating  from  a  vertex.  Consider  the 
difference  between  the  two  cross-boundary  derivatives  derived  from  two 
adjacent  microtriangles,  expressed  as  a  univariate  Bezier  polynomial. 


^  effected  by 
interpolation 
*  sect.  2.2. 1 
°  sect.  2.2.2 
°  sect.  2.2.3 


Figure  3:  First 
Stencils  a 


Order  Differentiability 
ad  their  Centers 


Interpolation  at  the  vertices  forces  the  value  of  that  difference  as  well  as 


its  tangential  derivative  to  vanish  at  the  vertex.  Thus  the  two  Bezier 
ordinates  closest  to  the  vertex  will  be  zero,  so  that  only  the  Bezier 
ordinates  shown  in  figure  3  are  affected.  The  total  number  of  conditions  so 
obtained  is  27. 

The  specific  stencil  illustrated  in  figure  3  means  that  algebraically 
3c00 10400  "  c0 110300  "  c0020300  ~  c0011300  “  0 

2.2.2.  C1  Conditions  on  Internal  Edges  from  Vertices  to  the  Centroid 

We  proceed  similarly  as  in  the  analysis  leading  to  the  equations  in  the 
preceding  subsection.  Consider  for  example  e14,  which  is  shared  by  the 
microtriangles  T14g  and  T147.  In  order  to  obtain  a  symmetric  stencil  we 
consider  the  derivative  in  the  direction  ®67  “  v7  ”  v6  which  we  express  on 
T^gg  as  6^^  m  ( 5eg4  t  ag^)/3  and  on  T^47  ns  ®67  *  ^®17  ^  Then 

we  differentiate  the  interpolant  on  each  microtriangle,  restrict  to  e14,  and 
consider  the  difference  a  univariate  Bezier  polynomial  in  b4  and  b ?.  The 
two  Bezier  ordinates  of  that  polynomial  that  are  closest  to  the  vertex  vanish 
because  of  interpolation  to  the  vertex  data,  and  the  one  at  the  centroid  will 


'j 

'wf 


be  dealt  with  in  the  next  subsection.  So  we  require  that  the  remaining  two 
vanish,  yielding  e.g.  the  condition 


10c 


2003000 


+  2c300200  “  6c2002010  ”  6c2002001 


This  is  the  particular  equation  depicted  in  figure  3.  On  the  other  edges 
we  proceed  analogously.  Altogether  we  obtain  six  conditions. 

2.2.3  C1  Condition  at  the  Centroid 

In  Alfeld,  1984,  it  was  shown  that  for  the  single  Clough-Tocher  split  the 


the  average  of  its  three  neighbors.  That  argument  does  not  apply  in  the 
present  case  because  the  interpolant  is  piecewise  polynomial  on  each  of  the 
subtriangles,  rather  than  polynomial.  However,  the  result  is  the  samel 

The  following  geometric  argument  (which  superseded  an  elaborate  algebraic 
one)  is  due  to  C.S.  Petersen  (private  communication):  For  the  interpolant  to 
be  C2  at  V4,  the  set  S,  say,  of  all  Bezier  control  points  corresponding 
to  cxxxixxx  where  i  >  3  and  the  x's  are  arbitrary  must  lie  in  the  same 
plane.  The  previously  imposed  conditions  imply  that  any  arrow  shaped  subset 
of  S  pointing  to  the  centroid  (like  {cooosooO'  c0 104000'  c00 14000 • 

c0004100*  or  ^c0005000'  c0 104000'  °0004100»  c0004001^  is  coplanar. 
Obviously,  all  points  in  S  will  be  coplanar  if  the  set  corresponding  to 

{c1004000 '  c0104000'  c0014000,  CqqosoOO^  (which  links  all  those  arrows)  is 
coplanar.  This  yields  the  condition 

3c0005000  “  C100400Q  "  c0 104000  "  c0014000  “  0 
as  depicted  in  figure  3. 

The  total  number  of  first  order  differentiability  conditions  is  33. 

2.3  Internal  Second  Order  Differentiability 

When  forcing  continuity  of  second  derivatives,  care  is  needed  that  the 
linear  system  describing  the  interpolant  does  not  become  overdetermined.  The 
mechanics  of  deriving  conditions  are  as  in  the  C1  case.  We  proceed  in 
several  stages. 

2.3.1  C2  Conditions  at  Vertices 

Second  order  differentiability  is  enforced  by  interpolation  to  the 
second  order  vertex  data. 
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2.3.2  C  Conditions  at  Subcsntroids 

Second  order  differentiability  is  serendipitously  implied  by  the 
first  order  conditions ,  for  details  see  Alfeld,  1984. 

2.3.3  C  Conditions  on  Lines  from  Vertices  to  Subcentroids 

Consider  the  difference  between  suitable  second  order  cross-boundary 
derivatives  as  they  have  been  obtained  from  two  neighboring  microtriangles. 
That  difference  is  a  univariate  cubic  polynomial  with  four  degrees  of 
freedom.  We  must  enforce  conditions  that  make  the  difference  zero.  One 
condition  each  has  been  inqtosed  at  the  vertex  and  at  the  subcentroid,  leaving 
two  degrees  of  freedom. 

For  example,  on  edge  we  consider  the  second  order  derivative  in  the 

direction  of  e14  and  obtain 

c3001001  “  e3 100001  “  c2200001  +  3c2100002  +  c2002001  “  3c2001002  “  0 
and  one  more  such  condition  that  can  be  obtained  by  shifting  the  stencil  of 

the  above  equation  towards  V7.  Altogether  there  are  12  conditions  of  this 
type. 

2 

2.3.4  C  Conditions  on  Lines  from  Subcentroids  to  the  Centroid 

Since  we  have  not  yet  enforced  a  second  order  condition  at  the  centroid, 
one  might  expect  that  three  conditions  are  needed.  However,  the  conditions 
described  in  section  2.3.3  imply  that  all  first  order  tangential  derivatives 
of  second  order  cross  boundary  derivatives  are  continuous  at  the  subcentroid. 

To  see  this  consider  e.g.  point  Vg  and  edge  e4g.  It  is  sufficient  to 
consider  any  particular  second  order  cross-boundary  derivative.  For 
convenience  we  choose  the  derivative 


Thus  we  are  interested  in  the  continuity  of  - - - - - -  at  V, . 

S*46a*168*36  * 

That  derivative  agrees  between  T34g  and  T136  since  it  is  a  tangential 
derivative  of  a  second  order  cross-boundary  derivative  in  the  direction  of 
•35  and  because  we  have  already  forced  second  order  differentiability 
across  e36.  Similarly  it  agrees  between  T14g  and  T147.  Hence  it  agrees 
between  *34$  and  T14g  across  «4g  which  is  what  we  wanted  to  show. 

Because  of  this  serendipitous  extra  piece  of  smoothness  we  have  to 
enforce  only  two  conditions  on  each  line  from  a  subcentroid  to  the  centroid. 
Considering  the  second  order  derivative  in  the  direction  of  e13  across 
e^4  we  obtain  the  symmetric  condition 

“^OOISOIO  “  c00220 10  +  3c00 12020  "  3e1002020  +  c2002010  +  c1003010  "  0 
This  particular  condition  is  depicted  in  figure  4,  and  there  are  six  such 

conditions  altogether. 

Note  that  the  conditions  derived  in  this  subsection  are  analogous  to 
those  in  the  previous  subsection,  which  is  as  one  would  expect  because  of  the 
symmetry  in  the  subtriangles.  However,  they  required  a  separate  derivation 
because  no  interpolation  takes  place  at  the  centroid. 

2.3.5  C2  Conditions  at  the  Centroid 

This  follows  serendipitous ly  as  at  the  subcentroids. 

2 

2.3.6  C  Conditions  on  Lines  from  Vertices  to  the  Centroid 

It  appears  at  first  that  as  in  section  2.3.3  two  conditions  have  to  be 
imposed  on  each  line,  yielding  six  conditions  altogether.  However,  one  of 
then  is  implied  by  the  other  five. 


Consider,  for  example,  the  second  order  e57  derivative  across  edge 
e24*  On  ^247  *•  can  write  e57  -  (e27  ♦  5047 )/3  and  on  T245  we  have 
057  “  (e52  +  5ej4>/3.  Carrying  out  the  nanipulations  yields  the  condition 

c0202100  ”  c020200 1  +  ^OIOSIOO  "  5c0103001  "  3c0102200  +  3c0102002  *  0 
which  is  depicted  in  figure  4.  The  number  of  conditions  so  obtained  equals 

six. 

Mow  consider  the  condition  with  cooi400  at  its  center.  Making  the 
corresponding  Besier  ordinate  of  the  difference  between  two  second  order 
derivatives  zero  is  equivalent  to  staking  the  tangential  derivative  at  the 
centroid  zero.  So  consider  for  example  the  derivative 


J*143*243e34  4 


and  assusw  that  we  have  enforced  the  other  five  Cr  conditions.  It  can  be 


seen  by  a  refinement  of  the  argustent  in  section  2.3.4  that  D  is  indeed 
continuous.  The  key  obeervation  is  that  the  relevant  triples  of  vertices, 
subcentroids,  and  centroid  are  colinear,  and  that  hence  D  is  a  tangential 
derivative  of  a  second  order  cross-boundary  derivative  in  each  microtriangle 
sharing  the  centroid. 

Thus  in  this  subsection  we  obtain  five  conditions  altogether. 


2.4  Intertriangular  Smoothness 

Here  we  must  consider  perpendicular  cross-boundary  derivatives.  The 
first  order  such  derivative  is  in  general  quartic  with  five  parameters,  but  we 
are  only  given  four  data.  So  we  require  additionally  that  the  leading 
coefficient  vanishes,  making  the  derivative  cubic.  Similarly,  we  require  that 


the  second  order  perpendicular  cross-boundary  derivative  be  linear 


A  smoothness  by 
interpolation 
£  smoothness  by 
serendipity 
♦  sect.  2.3.3 
°  sect.  2.3.4 
O  sect.  2.3.6 


e  that  of  o 


ier  Differentiability 
heir  Centers 
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It  is  not  necessary  to  consider  perpendicular  cross-boundary  directions 
of  unit  length.  Instead,  the  normalization  can  be  chosen  so  as  to  be 
appropriate  to  the  geometry  of  the  macrotriangle.  For  example,  on  edge  e^ 
we  define  the  normal  n12  by 

n12  “  *17  +  Y3*12  where  Y3  "  ”(e17*12)/(*17*12) 
and  express  the  cross-boundary  derivatives  as  univariate  polynomials  in  b^. 

Setting  the  leading  coefficient  of  the  first  order  cross-boundary  derivative 

equal  to  zero  yields  the  stencil  given  at  the  bottom  of  figure  5.  Similarly, 

setting  the  two  leading  coefficients  of  the  second  order  cross-boundary 

derivative  to  zero  yields  two  similar  conditions  which  are  depicted  in  figure 

5  on  other  edges.  All  twelve  such  conditions  can  be  obtained  from  figure  5  by 

rotating  the  labeling  and  modifying  appropriately  the  subscripts  of  Y . 

The  total  number  of  conditions  for  e  tratriangular  smoothness  is  twelve. 

2.5  Condensation  of  Parameters 

At  this  stage.  111  conditions  have  been  imposed,  leaving  10  degrees  of 
freedom.  In  disposing  of  the  ambiguity  one  would  like  to  meet  three 
objectives! 

1.  The  maximum  possible  decree  of  precision  (cubic)  is  maintained. 

2.  The  lnterpolant  is  independent  of  the  labeling  of  the  vertices. 

3.  The  additional  conditions  are  simple.  For  example,  one  would  not  like 
to  require  interpolation  to  additional  data  that  would  have  to  be  made  up. 

All  of  these  objectives  are  met  by  the  following  requirements: 

1.  The  lnterpolant  is  quartic  on  lines  from  vertices  to  the  centroid  (3 
conditions).  This  requirement  is  formulated  algebraically  by  setting  the 
fifth  order  tangential  derivatives  (which  are  constants)  equal  to  zero. 


2.  The  interpolant  is  cubic  on  lines  from  subcentroids  to  the  centroid  (6 
conditions).  Here  we  set  fourth  and  fifth  derivatives  to  zero  at  the 
centroid. 

3.  The  sum  of  tangential  fourth  order  derivatives  on  the  lines  from  the 
vertices  to  the  centroid  equals  zero.  More  precisely 


These  requirements  give  rise  to  particularly  simple  stencils  that  are 
depicted  in  figure  6. 

2.6  The  Solution  of  the  Linear  System 

The  analysis  in  the  previous  subsections  defines  a  linear  system  of  121 
equations  for  the  121  coefficients  of  the  interpolant.  That  system  was  set  up 
and  solved  using  the  symbol  manipulation  language  —MICE  (Hearn,  1983).  The 
first  48  equations  (corresponding  to  the  interpolation  conditions  and  the 
requirement  that  first  order  perpendicular  cross-boundary  derivatives  be 
quartic)  are  lower  triangular.  The  remaining  equations  were  reduced  to  lower 
triangular  form,  and  a  listing  of  the  coefficients  in  Forward  Elimination  form 
is  given  in  the  appendix. 

2.7  Precision  of  the  Interpolant. 

The  interpolant  developed  here  reproduces  all  bivariate  cubic  polynomials 
exactly.  This  follows  as  in  section  3.6  of  Alfeld,  1984. 


$  quartic 


o  sum  of  fourth 


Wgure  6: 


3.  Computational  Aspects. 

It  is  impractical  to  base  an  implementation  of  the  interpolant  on  the 
printed  formulas  given  in  the  appendix.  Instead,  any  code  should  be  written 
directly  by  the  symbol  manipulation  language  used  to  solve  the  linear 
system.  A  pilot  version  of  such  a  code  is  available  from  the  author.  It  code 
has  been  tested  by  the  analyzing  tool  MICROSOOPB  (Alfeld  and  Harris,  1984)  and 
it  was  verified  that  the  code  does  indeed  posses  the  smoothness, 
interpolation,  and  precision  properties  that  are  implied  by  the  mathematical 
construction. 

Conclusions 

The  scheme  developed  here  is  the  first  explicitly  given  piecewise 

2  2 
polynomial  C  interpolant  for  triangular  C  data. 
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Appendix:  Coefficients  of  Interpolant 


Following  are  expressions  for  the  coefficients  of  the  interpolant  as  the 
have  been  generated  by  HDUCS.  To  eliminate  typing  and  type-setting  errors 
the  results  have  been  reproduced  photographically.  The  organization  of  the 
listing  is  like  that  in  Alfeld,  1984.  The  parameters  have  been 

abbreviated  by  Gi. 
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C5000000  =  F{ 0  ,V1 ) 


2:  C4100000 

3:  C3200000 

4:  C4000001 

5:  C3000002 

6:  C3100001 

7s  c40 10000 
8:  C3020000 

9:  C4000010 

10:  C3000020 

11:  C3010010 

12:  c400 1000 

13:  C3002000 

14:  c300 10 10 

15:  C3001001 

16:  c0500000 

17:  C0410000 

18:  C0320000 

19:  c0400100 

20:  C0300200 

21:  C0310100 

22:  C1400000 

23:  C2300000 

24:  C04000Q1 

25:  C0300002 

26:  C1300001 

27: 


( 5*c5000000  +  F( 12 ,V1 ) )/5 

(  -  20*c5000000  +  40*c4100000  +  F( 12 12,V1 ) )/20 
( 5*c5000000  +  F(17,V1))/5 

(  -  20*c5000000  +  40*c400000 1  +  F( 1717, VI ) )/20 

(  -  20*c5000000  +  20*c4 100000  +  20*c4000001  +  F( 1217, VI ) )/20 

( 5*c5000000  +  F(13,V1))/5 

(  -  20*c5000000  +  40*c40 10000  +  F( 1313 , VI ) ) /20 
( 5*c5000000  +  F( 16, VI ) )/5 

(  -  20*c5000000  +  40*c40000 10  +  F( 1616, VI ) )/20 

(  -  20*c5000000  +  20*c4010000  +  20*c4000010  +  F( 13 16, VI ) )/20 

( 5*c5000000  +  F( 14, VI ) )/5 

(  -  20*c5000000  +  40*c400 1000  +  F( 1414,V1 ) )/20 

(  -  20*c5000000  +  20*c400 1000  +  20*c4000010  +  F( 1416,V1 ) )/20 

(  -  20*c5000000  +  20*c4001000  +  20*c4000001  +  F( 1417, VI ) )/20 

F(  0  ,V2 ) 

(5*c0500000  +  F(23,V2))/5 

(  -  20*c0500000  +  40*c0410000  +  F(2323,V2) )/20 
(5*c0500000  +  F(25,V2) )/5 

(  -  20*c0500000  +  40*c0400 100  +  F(2525,V2 ) )/20 

(  -  20*c0500000  +  20*c0410000  +  20*c0400100  +  F(2325,V2) )/20 

(5*c0500000  +  F( 2 1 ,V2 ) )/5 

(40*c1400000  -  20*c0500000  +  F( 2 12 1 , V2 ) )/20 
( 5*c0500000  +  F( 27, V2) )/5 

(  -  20*c0500000  +  40*c040000 1  +  F(2727,V2) )/20 
(20*c1400000  -  20*c0500000  +  20*c0400001  +  F(2127,V2) )/20 


C0401000  =  (5*c0500000  +  F(24,V2))/5 


28:  C0302000  -  (  -  20*c0500000  +  40*c0401000  +  F(2424,V2) )/20 

29:  C0301001  -  (  -  20*c0500000  +  20*c0401000  +  20*c0400001  +  F(2427,V2) )/20 

30:  C0301100  -  (  -  20*c0500000  +  20*c0401000  +  20*c0400100  +  F(2425,V2) )/20 

31:  cOOSOOOO  *  F<0,V3) 

32:  C1040000  -  (5*c0050000  +  F(31,V3))/5 

33:  C2030000  -  (40*c1040000  -  20*c0050000  +  F(3131,V3) )/20 

34:  C0040010  -  (5*c0050000  +  F(36,V3))/5 

35:  C0030020  -  (  -  20*c0050000  +  40*c0040010  +  F( 3636,V3) )/20 

36:  C1030010  -  (20*c1040000  -  20*c0050000  +  20*c0040010  +  F(3136,V3) )/20 

37:  C0140000  -  (5*c0050000  +  F(32,V3))/5 

38:  C0230000  -  (40*c0140000  -  20*c0050000  +  F(3232,V3) )/20 

39:  C0040100  -  (5*c0050000  +  F(35,V3))/5 

40:  C0030200  *  (  -  20*c0050000  +  40*c0040100  +  F( 3535, V3 ) )/20 

41:  cO 130100  -  (20*c0140000  -  20*c0050000  +  20*c0040100  +  F(3235,V3) )/20 

42:  C0041000  -  (5*c0050000  +  F<34,V3))/5 

43:  C0032000  -  (  -  20*c0050000  +  40*c0041000  +  F(3434,V3) )/20 

44:  C0031100  -  (  -  20*c0050000  +  20*c0041000  +  20*c0040100  +  F(3435rV3) )/20 

45:  c0031010  -  (  -  20*c0050000  +  20*c0041000  +  20*c0040010  +  F(3436,V3) )/20 

46:  C2200001  -  (c5000000*G3  +  oSOOOOOO  -  5*c4100000*G3  -  4*c4100000  - 

C4000001  +  10*c3200000*G3  +  6*c3200000  +  4*c3100001  -  10*c2300000*  G3  - 
4*c2300000  +  5*c1400000*G3  +  c1400000  +  4*c1300001  -  c0500000*  G3  -  c0400001) 
/6 

47:  C0220100  -  (c0500000*G1  +  C0500000  -  5*c0410000*G1  -  4*c0410000  - 

C0400100  +  10*c0320000*G1  +  6*c0320000  +  4*c0310100  -  10*c0230000*  G1  - 
4*c0230000  +  5*c0140000*G1  +  c0140000  +  4*c0130100  -  c0050000*  G1  -  c0040100) 
/6 

48:  C2020010  -  (  -  c5000000*G2  +  5*c4010000*G2  +  c4010000  -  C4000010  - 

10*c3020000*G2  -  4*c3020000  +  4*c3010010  +  10*c2030000*G2  +  6*  c2030000  - 
5*c1040000*G2  -  4*c1040000  +  4*c1030010  +  c0050000*G2  +  c0050000  -  c0040010)/6 
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49:  C0201002  -  ( 14661*c5000000*G3**2  +  29322*c5000000*G3  +  5535* 

c5000000*G2**2  +  16315*c5000000  -  22707*c4100000*G3**2  -  16092*  c4100000*G3  + 
6615*c4100000  -  16605*c40 10000*G2**2  -  11070*c4010000*  G2  -  5190*c4001000  + 
11070*c40000 10*G2  -  29322*c4000001*G3  -  29322*  c4000001  -  55782*c3200000*G3**2 

-  127656*c3200000*G3  -  57213*c3200000  +  16092*c3100001*G3  -  10665*c3 10000 1  + 

1 1070*c3020000*G2**2  +  22140*  c3020000*G2  +  5535*c3020000  -  22140*c3010010*G2  - 
13635*c3010010  +  5360*c3002000  +  4815*c3000020  +  15021*c3000002  + 

1 56978*c2300000*  G3**2  +  186300*c2300000*G3  +  35937*c2300000  + 
127656*c2200001*G3  +  98280*c2200001  +  11070*c2030000*G2**2  -  5535*c2030000  + 
5940*c20200 10  -  129087*c1400000*G3**2  -  7 1874*c1400000*G3  -  186300*c1300001*G3 

-  90585*c1300001  -  16605*c1040000*G2**2  -  22 140*c1040000*G2  -  5535*  c1040000  + 
22140*c1030010*G2  +  8505*c10300 10  +  35937*c0500000*G3**2  -  5643*c0500000*G1**2 

-  1 1286*c0500000*Gl  -  2021*c0500000  +  16443*  c0410000*G1**2  +  21600*c0410000*G1 
+  5157*c0410000  -  15030*c040 1000  +  11286*c0400100*G1  +  1 1286*c0400100  + 

7 1874*c040000 1 *G3  -  9342*  c0320000*Gl**2  +  2916*c0320000*G1  +  6615*c0320000  - 
21600*c0310 100*  G1  -  8235*c0310100  +  24392*c0302000  +  16632*c0301001  - 
5283*c0300200  +  36297*c0300002  -  14202*c0230000*G1**2  -  25488*c0230000*G1  - 

6129*  C0230000  -  2916*c0220100*G1  -  8586*c0220100  +  18873*c0140000*G1**2  + 

12258*c0140000*G1  +  25488*cOJ30100*G1  +  14823*c0 130 100  -  6129*  c0050000*G1**2  + 
5535*c0050000*G2**2  +  1 1070*c0050000*G2  +  7189*  c0050000  -  5l90*c0041000  - 
12258*c0040 100*G1  -  11070*c0040010*G2  -  11070*c0040010  +  5360*c0032000  - 
5769*c0030200  +  4815*c0030020 )/49896 

50:  C0004100  «  (  -  2681775*c5000000*G3**2  -  5363550*c5000000*G3  - 

810945*c5000000*G2**2  -  1594689*c5000000  +  4242915*c4100000*G3**2  + 

3 122280*c4 100000*G3  -  1120635*c4100000  +  2992005*c4010000*G2**2  + 

162 1890*c40 10000 *G2  -  3548930*c400 1000  -  1621890*c4000010*G2  + 

5363550*c400000 1*G3  +  5363550*c4000001  +  9846090*c3200000*G3**2  + 

228 14460*c3200000*G3  +  10286595*c3200000  -  3122280*c3100001*G3  + 

17 15445*c3 100001  -  3858570*c3020000*G2**2  -  4362120*c3020000*G2  - 
810945*c3020000  +  4362120*c3010010*G2  +  2147715*c3010010  +  3588200*  c3002000  - 
741825*c3000020  -  2834055*c3000002  -  281780 10*c2300000*G3**2  - 
33541 560*c2300000*G3  -  6484185*c2300000  -  22814460*c2200001*G3  - 
17822430*c2200001  +  1733130*c2030000*G2**2  +  3355020*c2030000*G2  + 

1 370 115*c20 30000  -  3355020*c20200 10*G2  -  2247750*c2020010  +  23254965* 
c1400000*G3**2  +  12968370*c1400000*G3  +  33541560*c1300001*G3  + 

16244955*c 1300001  +  196155*c1040000*G2**2  -  1 1 1240*c1040000*G2  - 
307395*c 1040000  +  111240*c1030010*G2  +  581445*c1030010  -  6484185* 
c0500000*G3**2  -  49815*c0500000*G1**2  -  99630*c0500000*G1  +  760707*  c0500000  + 
249075*c0410000*G1**2  +  398520*c0410000*Gl  +  149445*  c0410000  - 
21661 10*c040 1000  +  99630*c0400 100*G1  +  99630*c0400 100  -  12968370*c0400001*G3  - 
498150*c0320000*G1**2  -  597780*c0320000*Gl  -  149445*c0320000  - 
398520*c0310100*G1  -  165915*c0310 100  +  209840*  c0302000  -  3409560*c030 100 1  + 
33345*c0 300200  -  6636465*c0300002  +  498150  *c0230000*G1**2  +  398520*c0230000*Gl 
+  49815*c0230000  +  597780*  c0220100*Gl  +  4652 10*c0220100  +  10228680*c0201002  - 
249075*c0 140000*  Gl**2  -  99630*c0 140000*G1  -  398520*c0130100*Gl  - 
66285*c0 130100  +  49815*c0050000*G1**2  -  251775*c0050000*G2**2  - 
503550*c0050000*G2  +  962187*c0050000  -  4183310*c0041000  +  99630*c0040 100*G1  + 
503550*  c0040010*G2  +  503550*c0040010  +  41 1 1400*c0032000  +  132975*c0030200  - 
182655*c0030020  J/4532220 


-24- 


51 i  C0004010  -  <707670*c5000000*G3**2  +  1415340*c5000000*G3  ♦  154170 
*c5000000*G2**2  +  539058*c5000000  -  1197045*c4100000*G3**2  -  978750* 
c41 00000 *G3  +  2 18295*c4 100000  -  588195*c4010000*G2**2  -  308340*  c4010000*G2  + 
5436 10*c400 1000  ♦  308340*c4000010*G2  -  1415340*c4000001  *G3  -  1415340*c4000001 

-  2288520*c3200000*G3**2  -  5555790*c3200000*  G3  -  2559600*c3200000  + 

978750*c3 10000 1*G3  -  331425*c3100001  +  811080*  c3020000*G2**2  + 
868050*c3020000*G2  +  154170*c3020000  -  868050*  c3010010*G2  -  413505*c3010010  - 
553300*c3002000  +  142290*c3000020  +  740070*c3000002  +  6971 130*c2300000*G3**2  + 
8386470*c2300000*G3  +  1633635*c2300000  +  5555790*c2200001*G3  +  4403565*c2200001 

-  445770*  c2030000*G2**2  -  7541 10*c2030000*G2  -  279855*c2030000  +  754110* 
c2020010*G2  ♦  475065*c2020010  -  5826870*c1400000*G3**2  -  3267270*  c1400000*G3  - 
8386470* cl 30000 1 *G3  -  4088070*c1300001  +  40230*c1040000  *G2**2  + 
137430*c1040000*G2  +  97200*c1040000  -  137430*c1030010*G2  -  173880*c1030010  + 
1633635*c0500000*G3**2  -  103158*c0500000  +  216340*  c0401000  + 

3267270*c040000 1*G3  -  20520*c0310100  +  268820*c0302000  +  841320*c0301001  - 
20520*c0300200  +  1666035*c0300002  -  41040*c0220100  -  2523960*c0201002  - 
20520*c0130100  +  28485*c0050000*G2**2  +  56970*  c0050000*G2  -  168645*c0050000  + 
686200*c0041000  -  56970*c0040010*G2  -  56970*c0040010  -  670900*c0032000  - 

20520*c0 030200  +  16605*c0030020  +  1018710*c0004100 J/199260 

52 »  cOOl 1300  -  (  -  670491 00 *c5000000*G3**2  -  134098200*c5000000*G3  + 
155878560*c5000000*G2**2  -  1 1 10940473*c5000000  +  335245500*c4 100000*  G3**2  + 

536392800*c4100000*G3  +  20 1 147 300* c4 100000  -  497174220*  c4010000*G2**2  - 
311757120*c4010000*G2  +  345281 1905*c400 1000  ♦  31 1757120*c4000010*G2  + 

1 34098200*c400000 1 *G3  +  134098200*c4000001  -  67049 1000*c3200000*G3**2  - 
804589200*c3200000*G3  -  201 147300*c3200000  -  536392800*c3100001*G3  - 
320686560*c3 100001  +  42991 1280*c3020000*  G2**2  +  68259 1320*c3020000*G2  + 
155878560*c3020000  -  68259 1320*c30 100 10  *G2  -  393365160*c3010010  - 
3355729790*c3002000  +  136076220*c3000020  -  52490 160*c3000002  + 

67049 1000*c2300000*G3**2  +  536392800*c2300000*G3  +  67049100*c2300000  + 
804589200*c220000 1 *G3  +  431412480*c2200001  ♦  134525880*c2030000*G2**2  - 
177231240*c2030000*G2  -  185417100*c2030000  +  17723 1240*c2020010*G2  + 

237 1 56660*c20200 10  -  335245500*c1400000*  G3**2  -  134098200*c1400000*G3  - 
536392800*c1300001*G3  -  186588360*  C1300001  -  349481520*c1040000*G2**2  - 
446283000*c1040000*G2  -  96801480  *c1040000  +  446283000*c1030010*G2  + 
141533460*c1030010  +  67049100*  c0500000*G3**2  -  105026921 1*c0500000  + 

348470 1095*c040 1000  ♦  134098200*  c0400001*G3  +  22412970*c0310100  - 
3498749570*c0 302000  -  29656530*  C0301100  +  59733720*c0301001  + 

22412970*c0300200  +  81608040*c0300002  +  44825940*c0220100  -  179201160*c0201002 
+  224 12970*c0 130100  +  126340020*  c0050000*G2**2  +  252680040*c0050000*G2  - 
93556130 1*c0050000  +  3542861645*c0041000  -  252680040*c0040010*G2  - 
252680040*c00400 10  -  3496444550*c0032000  +  29656530*c0031 100  + 

224 12970*c00 30200  +  106537680*  c0030020  +  2658766140*c0004100  + 

824039730*c00040 10 )/5338 17540 


53:  c020 1200  -  ( 1288440*c5000000*G2**2  -  47455842*c5000000  +  4178520* 

c40 10000 *G2* *2  -  2576880*c40 10000*G2  +  1 57268345*c400 1000  +  2576880* 
c4000010*G2  -  1426 140*c3 10000 1  -  29598480*c3020000*G2**2  -  10933920* 
c3020000*G2  ♦  1288440*c3020000  +  10933920*c30 100 10*G2  -  1150740*  C3010010  - 
147384260*c3002000  +  462780*c3000020  -  1426140*c3000002  -  2852280*c220000 1  + 

508 39920 *c20 30000 *G2* *2  +  48263040*c2030000*G2  +  6755400*c2030000  - 
48263040*c20200 10*G2  -  18702360*c20200 10  -  1426140*  C1300001  - 
36040680*c1040000*G2**2  -  534 16800*c1040000*G2  -  17376120*  C1040000  + 
53416800*c1030010*G2  +  28134540*c1030010  -  48086379*c0500000  + 

16042 1030 *c040 1000  +  3852630*c03 10 100  -  164300 180*c0302000  +  3939840*  c0301100 

-  6366330*c030 100 1  +  3852630*c0300200  -  1426 140*c0300002  +  7705260*c0220 100  ♦ 
19098990*c020 1002  +  3852630*c0 130 100  +  9332280*  c0050000*G2**2  + 
18664560*c0050000*G2  -  37352694*c0050000  +  153414005*  c0041000  - 

1 8664560 *c00400 10 *G2  -  18664560*c00400 10  -  150400700*  C0032000  + 

2765070*c0031 100  +  3852630*c0030200  +  8506620*c0030020  -  49771260*c0011300  + 

1 60 374060* c0004 100  *  10 138095*c00040 10 )/20 1 14730 

54:  C1001003  -  (5749920*c5000000*G2**2  +  61904544*c5000000  -  21025440 

*c40 10000*G2**2  -  1 1499840*c40 10000*G2  -  207784690*c400 1000  +  11499840 
*c40000 10*G2  +  212 1255*c3 10000 1  +  26602560*c3020000*G2**2  +  30551040* 
c3020000*G2  +  5749920*c3020000  -  30551040*c3010010*G2  -  15275520*  c3010010  + 

2 10862 120*c3002000  +  1654425*c300 100 1  +  5149440*c3000020  +  2121255*c3000002  + 
42425 10*c220000 1  -  1 1 154240*c2030000*G2**2  -  22654080*c2030000*G2  - 
9525600*c2030000  +  22654080*c20200 10*G2  +  15275520*c20200 10  +  2121255*c1300001 

-  2147040*c1040000*G2**2  -  345600  *c1040000*G2  +  180 1440*c1040000  + 
345600*c10300 10*G2  -  3602880*  C1030010  +  60984738*c0500000  -  203185660*c0401000 

-  356535*c03 10 100  +  200364460*c0302000  -  3301 155*c030 1 100  +  1536435*c0301001  - 
356535*  C0300200  +  212 1255*c0300002  -  713070*c0220 100  +  8833860*c020 1200  - 
9572580*c020 1002  -  356535*c0 1 30 100  +  1974240*c0050000*G2**2  +  3948480* 
c0050000*G2  +  63516948*c0050000  -  205975510*c0041000  -  3948480*c00400 10  *G2  - 
3948480*c0040010  +  2 12278000*c0032000  +  356535*c0031 100  -  356535*  c0030200  + 
1373760*c00 30020  -  6417630*c001 1300  -  75277620*c0004100  -  88583490*c0004010) 
/29779650 

55:  C1010030  =  (9463392*c5000000*G2**2  -  45725472*c5000000  -  27897696 

*c40 10000*G2**2  -  18926784*c40 10000*G2  +  152992306*c400 1000  +  18926784 
*c4000010*G2  -  251883*c3 100001  +  16956864*c3020000*G2**2  +  36868608* 
c3020000*G2  +  9463392*c3020000  -  36868608*c30 100 10*G2  -  16152480*  C3010010  - 
1 548 18472*c300 2000  -  2281824*c3001010  -  240597*c300 100 1  +  10664352*c3000020  - 
251883*c3000002  -  503766*c220000 1  +  21881664*  c2030000*G2**2  + 
2954880*c2030000*G2  -  8970912*c2030000  -  2954880*  c2020010*G2  + 

22997952*c20200 10  -  251883*c1300001  -  30360096*c1040000*  G2**2  - 
40808448*c1040000*G2  -  10448352*c1040000  +  40808448*c1030010*  G2  + 
23178528*c10300 10  +  4330746*c100 1003  -  44303442*c0500000  +  145882156  *c0401000 
+  71307*c0310100  -  140412124*c0302000  +  405783*c030 1 100  -  225207*c030 100 1  + 

7 1307*c0 300200  -  25 1883*c0300002  +  142614*c0220100  -  1003428*c020 1200  + 

1 3974 12*c0 20 1002  +  71307*c0130100  +  9955872*  c0050000*G2**2  + 

1991 1744*c0050000*G2  -  35722404*c0050000  +  152756326*  c0041000  - 

1991 1744*c0040010*G2  -  1991 1744*c00400 10  -  155003152*  C0032000  -  71307*c0031 100 

-  2281824*c00310 10  +  71307*c0030200  +  11156832  *c0030020  +  1 283526*c00 1 1300  + 

9818820*c0004 100  +  133400898*c0004010)/41072832 


56:  C1002020  -  ( 12960*c5000000*G2**2  -  85392*c5000000  -  38880* 

c4010000*G2**2  -  25920*c4010000*G2  +  274190*c4001000  +  25920*  c4000010*G2  - 

2565*c3100001  +  25920*c3020000*G2**2  +  51840*c3020000*  G2  +  12960*c3020000  - 
51840*c3010010*G2  -  25920*c3010010  -  215480*  c3002000  +  2565*c3001001  + 
17280*c3000020  -  2565*c3000002  -  5130*  C2200001  +  25920*c2030000*G2**2  - 
1 2960 *c20 30000  +  25920*c20200 10  -  2565*c130000 1  -  38880*c1040000*G2**2  - 

51840*c1040000*G2  -  12960*  c1040000  +  51840*c1030010*G2  +  25920*c1030010  - 
46170*c1001003  -  76206  *c0500000  +  228260*c040 1000  +  2565*c03 10100  - 
150980*c0 302000  +  12825*  c0301100  -  12825*c0301001  +  2565*c0300200  - 
2565*c0300002  +  5130*  C0220100  -  30780*c020 1200  +  30780*c020 1002  + 

2565*c0130100  +  12960*  c0050000*G2**2  +  25920*c0050000*G2  -  84972*c0050000  + 

336890*  C0041000  -  25920*c0040010*G2  -  25920*c00400 10  -  395600*c0032000  -  2565 
*c0031 100  +  2565*c0030200  +  17280*c0030020  +  46170*c001 1300  -  1036260*  c0004100 
+  1676430*c0004010 )/492480 

57:  C0210200  -  (  -  285*c5000000  +  950*e4001000  -  950*c3002000  -  330* 

C0500000  +  1 175*c0401000  +  540*c0310100  -  1490*c0302000  -  486*c0301100  - 
54*c0301001  +  540*c0220 100  +  1458*c0201200  +  162*c0201002  -  285*  c0050000  + 
950*c004 1000  -  950*c0032000  +  675*c00040 10 )/1620 

58:  C0010040  -  (2790*c5000000*G2**2  -  16673*c5000000  -  6660*c4010000  *G2**2 

-  5580*c4010000*G2  +  55545*c400 1000  +  5580*c40000 10*G2  -  1260  *c3020000*G2**2  + 
7740*c3020000*G2  +  2790*c3020000  -  7740*c3010010*  G2  -  5010*c3010010  - 
55450*c3002000  -  570*c3001010  +  3 150*c3000020  +  15840*c2030000*G2**2  + 
10260*c2030000*G2  -  1080*c2030000  -  10260*  c2020010*G2  +  2160*c20200 10  - 
152 1 0*c1 040000 *G2**2  -  2 1420*c1040000*  G2  -  6210*c1040000  +  21420*c1030010*G2  + 
1 1850*c1030010  -  5130*  C1002020  -  16221*c0500000  +  53285*c0401000  - 
50930*c0302000  +  4500*  C0050000*G2**2  +  9000*c0050000*G2  -  12211*c0050000  + 
55735*c0041000  -  9000*c0040010*G2  -  9000*c0040010  -  55830*c0032000  + 
570*c0031010  +  6570*c00 30020  -  1710*c0004100  +  59760*c0004010)/15390 

59:  C0003200  -  (  -  138*c5000000  +  455*c4001000  -  440*c3002000  -  141* 

C0500000  +  470*c0401000  -  470*c0302000  -  141*c0050000  +  470*c0041000  - 
470*c0032000  +  675*c0004100)/270 

60:  C0010400  -  (  -  1620*c5000000*G2**2  +  62259*c5000000  -  25920* 

C4010000*G2**2  +  3240*c40 10000*G2  -  207530*c4001000  -  3240*c4000010*  G2  + 

1 19880*c3020000*G2**2  +  55080*c3020000*G2  -  1620*c3020000  -  55080*c3010010*G2  + 
3240*c3010010  +  207530*c3002000  -  2160*c3000020  -  187920*c2030000*G2**2  - 
184680*c2030000*G2  -  29160*c2030000  +  184680  *c2020010*G2  +  78840*c20200 10  + 
127980*c1 040000 *G2**2  +  191160*  c1040000*G2  +  63180*c1040000  - 
191 160*c1030010*G2  -  105840*c1030010  -  184680*c1010030  +  30780*c1002020  + 

60726*c0500000  -  199865*c040 1000  +  11970*c0310100  +  192200*c0302000  - 
1 1970*c0301 100  -  5130*c0300200  +  20520*c0220100  -  30780*c0210200  + 

6 1 560* c0 20 1200  +  8550*c0 130100  -  32400*c0050000*G2**2  -  64800*c0050000*G2  + 
26040*c0050000  *  188435*  C0041000  +  64800*c00400 10*G2  +  64800*c00400 10  + 

162500 *c0032000  -  8550  *c0031100  -  20520*c0031010  +  15390*c0030200  - 
63720*c0030020  +  277020*  c0011300  +  277020*c0010040  +  530955*c0004100  - 
361530*c0004010  -  548910*  c0003200 )/184680 


61:  C1001030  -  (37980*c5000000*G2**2  -  240573*c5000000  -  87840* 

c4010000*G2**2  -  75960 *c40 100 00 *G2  +  801910*c4001000  +  75960*c40000 10  *G2  - 
28440*c3020000*G2**2  +  99720*c3020000*G2  +  37980*c3020000  -  99720*c30 100 10*G2  - 
75960*c3010010  -  801910*c3002000  +  43920*c3000020  +  232560*c2030000*G2**2  + 

1 56600 *c20 300 00 *G2  -  1 1880*c2030000  -  156600*c20200 10*G2  +  9720*c20200 10  - 
2 18340*c1 040000 *G2**2  -  308520*  c1040000*G2  -  90 180*c1040000  + 
308520*c1030010*G2  +  166320*c1030010  +  126360*c1010030  -  5940*c1002020  - 

234042*c0500000  +  769255*c040 1000  -  8190*c0310100  -  736600*c0302000  + 

8190*c0301 100  +  3510*c0300200  -  14040  *c0220100  +  21060*c0210200  - 
42 120*c020 1200  -  5850*c0130100  +  64080*  c0050000*G2**2  +  128160*c0050000*G2  - 
173880*c00 50000  +  788845*  c0041000  -  128 160*c00400 10*G2  -  128 160*c00400 10  - 
771 100*c00 32000  +  5850*c0031 100  +  14040*c00310 10  -  1 0530*c0030200  + 

1 06200 *c0 030 020  -  189540*c001 1300  +  126360*c0010400  -  325620*c0010040  - 
363285*c0004100  +  915030*c0004010  +  375570*c0003200 J/90720 

62:  cl  100003  *  (23760*c5000000*G2**2  -  1 168 17*c5000000  -  47520* 
c40 10000 *G2**2  -  47520*c40 10000*G2  +  398680*c400 1000  +  47520*c40000 10  *G2  + 
15930*c3100001  -  47520*c3020000*G2**2  +  47520*c3020000*G2  +  23760*c3020000  - 
47520*c30 100 10*G2  -  23760*c3010010  -  434470*  C3002000  -  23760*c300 10 10  - 

15930*c300 100 1  +  23760*c3000020  -  270*  c3000002  +  31860*c220000 1  + 
190080*c2030000*G2**2  +  142560*c2030000*  G2  -  142560*c2020010*G2  + 

23760*c20200 10  +  15930*c1300001  -  166320*  c1040000*G2**2  -  237600*c1040000*G2  - 
7 1280*c 1040000  +  237600*  c1030010*G2  +  142560*c10300 10  -  213840*c10 10030  - 

7 1280*c 1002020  +  213840*c100 1030  +  140940*c100 1003  -  1 12332*c0500000  + 

376255*c040 1000  +  8100*c0310100  -  384400*c0302000  -  8100*c0301100  - 
15930*c030 1001  -  270*c0300002  +  8100*c0220 100  -  24300*c02 10200  +  24300*c020 1200 
+  46980  *c020 1002  +  47520*c0050000*G2**2  +  95040*c0050000*G2  -  65076*  c0050000 
+  377575*c0041000  -  95040*c00400 10*G2  -  95040*c0040010  -  384340*c0032000  + 
47520*c0030020  -  38205*c0004100  +  203580*c0004010  +  12 1770*c0003200 )/145800 

63:  C1000040  =  (95040*c5000000*G2**2  +  1 16817*c5000000  -  332640* 

c4010000*G2**2  -  190080*c4010000*G2  -  398680*c400 1000  +  190080*  c4000010*G2  - 
15930*c3 100001  +  380 160*c3020000*G2**2  +  475200*  c3020000*G2  +  95040*c3020000  - 
475200*c ( 30 100 10*G2  -  213840*c3010010  +  434470*c3002000  +  23760*c3001010  + 
15930*c300 100 1  +  166320*c3000020  +  270*c3000002  -  31860*c220000 1  - 
95040*c2030000*G2**2  -  285120*  c2030000*G2  -  142560*c2030000  + 

285120*c2020010*G2  +  261360*c2020010  -  15930*c130000 1  +  145800*c1 100003  - 
95040*c1 040000 *G2**2  -  95040*  c1040000*G2  +  95040*c10300 10*G2  +  213840*c1010030 
+  213840*c1001030  -  140940*c100 1003  +  1 12332*c0500000  -  376255*c0401000  - 
8100*c0310100  +  384400*c0302000  +  8100*c030 1 100  +  15930*c0301001  +  270*c0300002 
-  8100  *c0220 100  +  24300*c02 10200  -  24300*c020 1200  -  46980*c020 1002  +  47520* 
c0050000*G2**2  +  95040*c0050000*G2  +  160 1 16*c0050000  -  377575*  c0041000  - 
95040*c00400 10*G2  -  95040*c00400 10  +  384340*c0032000  +  47520  *c0030020  +  38205 
*c0004100  -  203580*c00040 10  -  12 1770*c0003200 J/641520 

64:  C0101300  *  (34065*c5000000  -  1 13550*c400 1000  +  1 13550*c3002000  + 

34054*c05 00000  -  1 13495*c040 1000  -  330*c0310100  +  1 13440*c0302000  +  330 
*c030 1 100  +  2970*c0300200  +  5940*c0210200  +  330*c0130100  +  34076*  c0050000  - 
1 13605*c0041000  +  1 1 3660*c0032000  -  330*c0031100  -  990*  C0030200  - 
66420*c00 11300  +  35640*c0010400  -  334035*c0004100  +  990*  c0004010  + 

3023 10*c0003200 )/48600 


65 i  cO 103100  -  (45*c5000000  -  150*c4001000  +  150*c3002000  ♦  46*  cOSOOOOO  - 
155*c0401000  +  160*c0 302000  +  44*c0050000  -  145*c0041000  ♦  140*c0032000  - 
75*c0004100  -  90*c0004010  +  90*c0003200 )/60 

66:  cO 102200  -  (261*c5000000  -  870*c4001000  +  870*c3002000  +  261*  cOSOOOOO  - 
870* c040 1000  +  870*c0302000  +  22*c0300200  +  33*c0210200  -  333*c0101300  + 
261*c0050000  -  870*c004 1000  +  870*c0032000  -  11*  C0030200  -  531*c0011300  + 

297* cOO 10400  -  2538*c0004100  ♦  231 1*c0003200)/33 

67:  C0022100  -  ( 1 161*c5000000  -  3870*c4001000  +  3870*c3002000  ♦  1161* 

C0500000  -  3870*c0401000  +  3870*c0302000  +  99*c0300200  +  11*c0220100  + 

165*c02 10200  +  33*c0201200  +  11*c0130100  -  66*c0102200  -  1611*  c0101300  + 

1 161*c0050000  -  3870*c0041000  +  3870*c0032000  -  11*c0031100  -  33*c0030200  - 
2205*c001 1300  +  1188*c0010400  -  11259*c0004100  +  10206  *c0003200 )/1 1 

68:  cO 120200  -  (  -  129*c5000000  +  430*c4001000  -  430*c3002000  -  129* 

C0500000  +  430*c040 1000  -  430*c0302000  -  11*c0300200  -  22*c0210200  +  11 
*cO 102200  +  168*c0101300  -  129*c0050000  +  430*c0041000  -  430*c0032000  + 

234*c001 1300  -  99*cOO 10400  +  1251*c0004100  -  1 134*c0003200 )/1 1 

69:  C2002010  -  (88851*c5000000  -  300815*c4001000  -  7965*c3100001  + 

318710*c3002000  +  7965*c3001001  +  135*c3000002  -  15930*c2200001  -  7965 
*c1300001  +  72900*c1 100003  -  70470*c1001003  +  87285*c0500000  -  292985*  c0401000 
-  4050*c0310100  +  300440*c0 302000  +  4050*c0301 100  +  7965*  c0301001  + 
135*c0300002  -  4050*c0220 100  +  12150*c0210200  -  12150*  c0201200  - 
23490*c0201002  -  40590*c0103100  +  86064*c0050000  -  286880*  c0041000  + 
286880*c0032000  -  31635*c0004100  -  162675*c0004010)/1 1880 

70:  c2 100002  -  (  -  77859*c5000000  +  2625 10*c400 1000  +  3510*c3100001  -  273340 
*c3002000  -  3510*c3001001  -  540*c3000002  +  5670*c2200001  +  5670  *c2002010  + 
2160*c1300001  -  24300*c1 100003  +  38880*c1001003  -  76770*  c0500000  + 

257065*c040 1000  -  260560*c0302000  -  2160*c0301001  -  540*  c0300002  + 
8910*c0201002  +  46260*c0103100  -  76071*c0050000  ♦  253570*  c0041000  - 
253570*c0032000  +  2340*c0004100  +  166725*c0004010)/4050 

71:  C0003020  ■  (  -  276*c5000000  +  920*c4001000  -  920*c3002000  -  276* 

C0500000  +  920*c0401000  -  920*c0302000  +  180*c0103100  -  273*c0050000  + 
905*c0041000  -  890*c0032000  +  225*c0004100  +  945*c0004010  -  270*  c0003200)/270 

72:  C1002002  -  (3093*c5000000  -  9106*c4001000  +  1566*c3100001  +  4972* 

C3002000  -  1566*c3001001  -  540*c3000002  +  1566*c2200001  -  4698*c2100002  + 
1566*c2002010  -  4860*c1 100003  +  23328*c1001003  +  3354*c0500000  -  1041 1*c0401000 
+  8104*c0302000  -  540*c0300002  +  3078*c0201002  +  3354*  C0050000  - 
1041 1*c004 1000  +  8104*c0032000  -  45657*c0004100  -  49572*  c0004010  + 
41526*c0003200  +  41526*c0003020 )/7776 

73:  cO 10 1003  -  (  -  573*c5000000  +  1574*c4001000  -  54*c3 100001  -  548* 

C3002000  +  54*c3001001  +  162*c3000002  -  54*c2200001  +  162*c2 100002  -  54 
*c2002010  +  1458*c1 100003  -  162*c1002002  -  2916*c1001003  -  582*  c0500000  + 
1619*c0401000  -  656*c0302000  +  162*c0300002  +  324*c0201002  -  582*c0050000  + 
1619*c0041000  -  656*c0032000  +  18819*c0004100  +  18954*  C0004010  - 
17334*c0003200  -  17334*c0003020)/3402 


74s  C0202100  -  ( 1860*c5000000  -  6329*c400 1000  +  67 16*c3002000  -  594* 

C2200001  -  1782*c2 100002  -  594*c1300001  +  5346*c1 100003  -  1782*c1002002  + 
486*c1001003  +  1959*c0500000  -  6824*c040 1000  +  7904*c0302000  +  594*  c0301001  - 
4860*c0 101003  +  1860*c0050000  -  6329*c004 1000  +  6716*  C0032000  +  5832*c0004 100 
+  4347*c00040 10  -  6966*c0003200  -  6966*  c0003020 )/594 

75s  c200 1002  *  ( 8889*c5000000  -  26228*c400 1000  +  16022*c3002000  -  1782 

*c3000002  -  594*c220000 1  -  1782*c2 100002  -  594*c1300001  -  10692*  cl  100003  + 
32562*c 100 1003  +  8988*c0500000  -  26723*c0401000  +  17210*  C0302000  + 

594*c0301001  -  1782*c0300002  -  594*c0202100  -  3564*c020 1002  +  32562*c0101003  + 
8889*c00 50000  -  26228*c0041000  +  16022*c0032000  -  202662*c0004 100  - 
204 147* cO 0040 10  +  183708*c0003200  +  183708*c0003020)/1782 

76s  cOI 10300  =  (  -  129*c5000000  +  430*c4001000  -  430*c3002000  -  129* 

C0500000  +  430*c040 1000  -  430*c0302000  -  11*c0300200  -  11*c0210200  -  11 
*c020 1200  +  201*c0101300  -  129*c0050000  +  430*c0041000  -  430*c0032000  + 

234*c001 1300  -  99*c0010400  +  1251*c0004100  -  1 134*c0003200 )/33 

77s  C0002030  =  (33*c5000000  -  110*c4001000  +  110*c3002000  +  9*  c1000040  +  33 

*c0500000  -  1 10*c040 1000  +  110*c0302000  +  33*c0050000  -  110*c0041000  + 

1 10*c0032000  +  9*c00 10040  -  369*c0004010  +  486*  c0003020 )/234 

78s  cO 100004  -  (51*c5000000  -  149*c4001000  +  86*c3002000  -  11*  C3000002  -  11 
*c2100002  -  1 1*c2001002  -  33*c1100003  +  201*c1001003  +  51  *c0500000  - 

149*c040 1000  +  86*c0302000  +  234*c0101003  +  51*c0050000  -  149*c0041000  + 
86*c0032000  -  1251*c0004100  -  1251*c00040 10  +  1134*  c0003200  +  1 134*c0003020 ) 
/99 

79s  C0002300  *  (3*c5000000  -  10*c4001000  +  10*c3002000  +  3*c0500000  - 

10*c0401000  +  10*c0302000  +  3*c0101300  +  3*c0050000  -  10*c0041000  +  10*c0032000 
+  3*c001 1300  -  36*c0004100  +  54*c0003200)/33 

80s  cOO 1 1030  *  (  -  3*c5000000  +  10*c4001000  -  10*c3002000  -  3*  c1001030  -  3 
*c0500000  +  10*c040 1000  -  10*c0302000  -  3*c0050000  +  10*  c0041000  -  10*c0032000 
+  36*c0004010  -  54*c0003020  +  33*c0002030 )/3 

81s  C0002003  =  ( 18*c5000000  -  61*c4001000  +  64*c3002000  +  3*c1001003  + 

18*c0500000  -  61*c040 1000  +  64*c0302000  +  3*c0101003  +  18*c0050000  - 

61*c0041000  +  64*c0032000  +  36*c0004100  +  36*c0004010  -  54*c0003200  - 

54*c0003020)/33 

82s  C0012200  =  -  c0102200  -  c0003200  +  3*c0002300 

83s  cOO 12020  =  -  c1002020  -  c0003020  +  3*c0002030 

84s  cO 102002  =  (  -  17*c5000000  +  57*c4001000  -  58*c3002000  -  18*  c1002002  - 
17*c0500000  +  57*c040 1000  -  58*c0302000  -  17*c0050000  +  57*  c0041000  - 

58*c0032000  +  18*c0003200  +  18*c0003020  +  54*c0002003 )/18 

85s  C0013100  =  -  C0103100  -  C0004100  +  3*c0003200 

86s  C1003010  =  ( 276*c5000000  -  920*c4001000  +  920*c3002000  +  276*  c0500000  - 

920*c040 1000  +  920*c0302000  +  273*c0050000  -  905*c0041000  +  890*c0032000  + 
90*c0013100  -  495*c0004100  -  585*c0004010  +  270*  c0003020)/90 


87 i  C0003002  -  (531*c5000000  -  1775*c4001000  +  1790*c3002000  -  90*  c1003010 
+  531*c0500000  -  1775*c0401000  +  1790*c0302000  -  90*c0103100  +  528*c0050000  - 
1760*c0041000  +  1760*c0032000  -  585*c0004100  -  585*  c0004010)/270 

88 i  C0004001  -  (3*c5000000  -  10*c4001000  +  10*c3002000  +  3*c0500000  - 
10*e0401000  +  10*c0302000  +  3*c0050000  -  10*c0041000  +  10*c0032000  -  3*c0004100 
-  3*c0004010)/3 

89:  C1004000  -  (  -  15*c5000000  +  50*c4001000  -  50*c3002000  -  15*  c0500000  + 

50*c0401000  -  50*c0 302000  -  15*c0050000  +  50*c0041000  -  50*  C0032000  + 

27*c00040 10  +  27*c0004001)/9 

90:  cO 104000  -  {  -  3*c5000000  +  10*c4001000  -  10*c3002000  -  9*  c1004000  -  3 
*c0500000  +  10*c0401000  -  10*c0302000  -  3*c0050000  +  10*  c0041000  -  10*c0032000 
♦  27*c0004001 )/9 

91:  C0005000  -  <3*c5000000  -  10*c4001000  +  10*c3002000  +  3*c0500000  - 

10*c040 1000  +  10*c0302000  +  3*c0050000  -  10*c0041000  +  10*c0032000)/9 

92:  C1003001  -  (c5000000  -  5*c4001000  +  10*c3002000  +  55*c1004000  -  30 

*c1003010  -  C0005000 )/30 

93:  cO 103001  -  (c0500000  -  5*c0401000  +  10*c0302000  +  55*c0 104000  -  30 
*c0103100  -  c0005000 )/30 

94:  C0013010  -  (33*c5000000  -  110*c4001000  +  110*c3002000  -  55*  c1004000  + 

33*c0500000  -  1 10*c0401000  +  110*c0302000  -  55*c0 104000  +  34  *c0050000  - 

1 15*c0041000  4  120*c0032000  -  30*c0013100  +  65*c0005000)/30 


95:  C2002001  -  (c5000000  -  5*c4001000  +  12*c3002000  -  6*c2002010  +  5* 

C1004000  -  c0005000)/6 

96:  C0202001  »  (c0500000  -  5*c0401000  +  12*c0302000  -  6*c0202100  +  5* 

cO 104000  -  c0005000)/6 

97:  C0022010  -  (3*c5000000  -  10*c4001000  4  10*c3002000  -  5*c1004000  + 

3*c0500000  -  10*c0401000  +  10*c0302000  -  5*c0104000  +  4*c0050000  -  15*c0041000 
4-  22*c0032000  -  6*c0022100  +  5*c0005000)/6 

98:  cO 100400  -  (c0300200  +  C0210200  +  c0201200  +  3*c0110300  +  3*  c0101300) 
/9 

99:  C0021200  -  -  c0120200  -  3*c0110300  -  c0030200  -  3*c0011300  +  9* 

cOO 10400 

100:  cO 200300  -  (c0300200  +  c0210200  4-  c0201200)/3 
101:  C0020300  -  (c0120200  +  c0030200  4  c0021200)/3 


102s  C2001020  *  (  -  2*c5000000*G2**2  +  7*c4010000*G2**2  +  4*  c4010000*G2  - 

4*c40000 10*G2  -  8*c3020000*G2**2  -  10*c3020000*G2  -  2*c3020000  +  10*c30 100 10*G2 
+  4*c30 100 10  -  5*c3000020  +  2*c2030000  *G2**2  +  6*c2030000*G2  +  3*c2030000  - 
6*c20200 10*G2  -  6*c2020010  +  2*c1040000*G2**2  +  2*c1040000*G2  -  2*c1030010*G2  - 
9*c10 10030  -  9  *c100 1030  +  27*c1000040  -  c0050000*G2**2  -  2*c0050000*G2  - 
cOOSOOOO  +  2*c00400 10*G2  +  2*c0040010  -  c0030020)/3 

103s  C0021020  =  (  -  c5000000*G2**2  +  2*c40 10000*G2**2  +  2*c4010000*  G2  -  2 

*c4000010*G2  +  2*c3020000*G2**2  -  2*c3020000*G2  -  c3020000  +  2*c30100 10*G2  + 
2*c3010010  -  C3000020  -  8*c2030000*G2**2  -  6*  c2030000*G2  +  6*c2020010*G2  + 
7*c1040000*G2**2  +  10*c1040000*G2  +  3*c1040000  -  10*c10300 10*G2  -  6*c1030010  - 
9*c1010030  -  2*c0050000  *G2**2  -  4*c0050000*G2  -  2*c0050000  +  4*c00400 10*G2  + 
4*c00400 10  -  5*c0030020  -  9*c0011030  +  27*c00 10040 )/3 

104s  C2000030  =  (2*c5000000*G2**2  -  7*c40 10000*G2**2  -  4*c4010000*  G2  +  4 

*c4000010*G2  +  8*c30 20000 *G2**2  +  10*c3020000*G2  +  2*  c3020000  -  10*c30 100 10*G2 
-  4*c3010010  +  5*c3000020  -  2*c2030000*  G2**2  -  6*c2030000*G2  -  3*c2030000  + 
6*c2020010*G2  +  6*c2020010  +  3  *c2001020  -  2*c1040000*G2**2  -  2*c1040000*G2  + 
2*c1030010*G2  +  c0050000*G2**2  +  2*c0050000*G2  +  c0050000  -  2*c00400 10*G2  -  2* 
C0040010  +  c0030020 )/9 

105s  C0020030  =  (c5000000*G2**2  -  2*c40 10000*G2**2  -  2*c40 10000*G2  + 

2*c40000 10*G2  -  2*c3020000*G2**2  +  2*c3020000*G2  +  C3020000  -  2*  c3010010*G2  - 
2*c3010010  +  C3000020  +  8*c2030000*G2**2  +  6*c2030000*  G2  -  6*c20200 10*G2  - 
7*c1 040000 *G2* *2  -  10*c1040000*G2  -  3*  c1040000  +  10*c10300 10*G2  +  6*c1030010  + 
2*c0050000*G2**2  +  4*  c0050000*G2  +  2*c0050000  -  4*c00400 10*G2  -  4*c0040010  + 
5*c0030020  +  3*c0021020 )/9 

106s  C1000004  -  (C3000002  +  c2100002  +  C2001002  +  3*c1 100003  +  3*  C1001003) 

/9 

107s  C1200002  -  -  3*c1 100003  -  c0300002  -  c0201002  -  3*c0101003  +  9* 

cO 100004 

108$  c2000003  =  (C3000002  +  c2100002  +  c2001002)/3 

109$  C0200003  =  (C1200002  +  c0300002  +  c0201002)/3 

110$  c0014000  »  (3*c5000000  -  10*c4001000  +  10*c3002000  -  5*c1004000  + 

3*c0500000  -  10*c040 1000  +  10*c0302000  -  5*c0104000  +  3*c0050000  -  10*c0041000 
+  10*c0032000  +  6*c0005000)/5 

Ills  cOOO 1040  =  -  C0005000  +  4*c0004010  -  6*c0003020  +  4*c0002030 

112s  cOOO 1400  =  -  cOOOSOOO  +  4*c0004100  -  6*c0003200  +  4*c0002300 

113s  C0000005  =  -  4*c0005000  +  15*c0004001  -  20*c0003002  +  10*  c0002003 

114s  C2003000  =  ( cSOOOOOO  -  5*c4001000  +  10*c3002000  +  5*c1004000  -  c0005000) 

/10 

115s  C0000050  =  C0005000  -  5*c0004010  +  10*c0003020  -  10*c0002030  +  5* 

cOOO 1040 
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116s  C0023000  -  (cOOSOOOO  -  5*c0041000  +  I0*c0032000  +  5*c0014000  -  c0005000) 

/10 

117*  cOOOOSOO  -  cOOOSOOO  -  5*c0004100  +  10*c0003200  -  10*c0002300  +  5* 
cOO0 1400 


11®*  cO 203000  -  (c0500000  -  5*c0401000  +  10*c0302000  +  5*c0104000  -  c0005000) 
/10 

119t  C0001004  -  (  -  cOOOSOOO  +  5*c0004001  -  10*c0003002 •  +  10*c0002003  + 

C0000005 )/5 

120j  C1020020  -  (c5000000*G2**2  -  2*c4010000*G2**2  -  2*c4010000*G2  + 

2*c40000 10*G2  -  2*c3020000*G2**2  +  2*c3020000*G2  +  c3020000  -  2*  c3010010*G2  - 
2*c3010010  +  C3000020  +  8*c2030000*G2**2  +  6*c2030000*  G2  -  6*c2020010*G2  - 
7*c1 040000 *G2**2  -  10*c1040000*G2  -  3*  C1040000  +  10*c1030010*G2  +  6*c1030010  + 
2*c0050000*G2**2  +  4*  c0050000*G2  +  2*c0050000  -  4*c0040010*G2  -  4*c0040010  + 
2*c0030020)/3 

121 i  c20 10020  -  (c5000000*G2**2  -  5*c4010000*G2**2  -  2*c4010000*G2  + 
2*c4000010*G2  +  10*c3020000*G2**2  +  8*c3020000*G2  +  c302000Q  -  8  *c3010010*G2  - 
2*c3010010  +  C3000020  -  10*c2030000*G2**2  -  12*  c2030000*G2  -  3*c2030000  + 
12*c2020010*G2  +  6*c2020010  +  5*c1040000  *G2**2  +  8*c1040000*G2  +  3*c1040000  - 
8*c1030010*G2  -  6*c1030010  +  3*c1020020  -  c0050000*G2**2  -  2*c0050000*G2  - 
C0050000  +  2*c0040010*  G2  +  2*c0040010  -  c0030020)/3 
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